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^— ^ i Abstract. It is proved that if the scattering amplitudes for two obstacles (from a large class of obstacles) differ a 

^—•^ , little, then the obstacles differ a little, and the rate of convergence is given. An analytical formula for calculating 

^v^j ' the characteristic function of the obstacle is obtained, given the scattering amplitude at a fixed frequency. 

a 

f — Introduction. 

Let D C t 3 be a bounded domain with a smooth boundary T, 

(V 2 + k 2 )u = in D' := R 3 \ D, k = const > 0; u = on T (1) 

u = exp(ika ■ x) + A(a', a, fc)r _1 exp(ikr) + o(r _1 ), r := \x\ — > oo, a' := xr^ 1 . (2) 

Here a £ S 2 is a given unit vector, S 2 is the unit sphere in R 3 , the function A(a',a,k) is called the 
scattering amplitude (the radiation pattern). It is well known [1] that problem (l)-(2) has a unique 
solution, the scattering solution, so that the map Y — * A(a' , a, k) is well defined. We consider the inverse 
obstacle scattering problem (IOSP): given A(a! , a) := A{a' 1 a, k = 1) for all a' , a £ S 2 and a fixed k (for 
example, take k = 1 without loss of generality), find Y. 
. Let us assume that Y C 7a, where 7> is the set of star-shaped (with respect to a common point O) 

surfaces, which are located in the annulus < ao < \x\ < a±, and whose equations X3 = (j)(xx,%2) in the 
local coordinates (in which X3 is directed along the normal to T at a point seT), have the property 

H\\c.»<co, (3) 
C 2 ' A is the space of twice diffcrcntiablc functions, whose second derivatives satisfy the Holder condition of 
order < A < 1, A and Co are independent of <p and Y. 

Uniqueness of the solution to IOSP with fixed frequency data is first proved in [1, p. 85]. We are 
■ interested here in the stability problem: suppose that Yj £ generate Aj{a! a), j = 1, 2, and 

max \Ax(a',a) - A 2 (a',a)\ < S. (4) 

a' ,a£S 2 

What can one say about the Hausdorff distance between D\ and D2'. p '■— sup 2 , e p 1 inf yg r 2 \x — y\. Let D\ 
denote a connected component of D\ \ D 2 , -D12 := D\ U D 2 , Ti2 := dD 12 , D' 12 := ^ 3 \ ^12, Ti := dD\ := 
Y[ U Y 2 , Y 2 C Y 2 := dD 2 , Y[ C T\ := dD\. Let us assume, without loss of generality, that p — \xq — j/o|> 
xq £ Y[, yo £ Y 2 . Can one obtain a formula for calculating Y, given A(a',a) for all a', a £ S 2 , k = 1 is 
fixed? No such formula is known for IOSP. For inverse potential scattering problem with fixed-energy data 
such a formula and stability estimates are obtained in [2] , [3] . These results are based on the works [7] , [8] , 
[10]-[17], [19]-[21]. 

In section II we prove that p < c\ f ^jfo^ J as S — * 0. We also prove some inversion formula, but it 

is an open problem to make an algorithm out of this formula. In Remark 3, we comment on some recent 
papers [4-6] in which attempts are made to study the stability problem and point out a number of errors 
in these papers. Our result, formulated as Theorem 1 in section II, is stronger than the results announced 
in Theorem 1 in [4], Theorem 1 in [5] and Theorem 2.10 in [6]. 
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II. Stability Result and a Reconstruction Formula. 

( In I In <5I \ C2 

Theorem 1. Under the assumptions of section I, one has p(S) < C\ ( "| n " 1 j , where c\ and c 2 are 
positive constants independent of S. 



Proposition 1. There exists a function v e (a,9) e L 2 (S 2 ) such that 

-4tt lim f A(6', a)u e (a, 6)da = ~xd(X). (5) 

e^0J S 2 2 

f l, x e D 

Here A e M 3 is an arbitrary fixed vector, Xd(x) := < , XdM '■= L 3 exp(— i\ ■ x)xD{x)dx, 

y o, x f. d 

9, 9' e M := {9 : 9 e C 3 , 9 -9=1}, 0' — 6 = X, and A(9' , a) is defined by the absolutely convergent series 

oo 

A{6',a) = Y j A e {a)Y l {6'), 9' e M, A e (a) := A(a> , a)Y^)da' , (6) 
to Js2 

where Yi(a) are the orthonormal in L 2 {S 2 ) spherical harmonics, Yg(6') is the natural analytic continuation 
of Ye(a') from S 2 to M, and the series (6) converges absolutely and uniformly on compact subsets of 
S 2 xM. 



Remark 1. The stability result given in Theorem 1 is similar to the one in [3], p. 9, formula (2.42), for 
inverse potential scattering. 

Remark 2. Proposition 1 claims the existence of the inversion formula (5). An open problem is to construct 
the function v e (a, 9) algorithmically, given the data A(a', a) Va', a e S 2 . 

Proof of Theorem 1. First, we prove that p(6) — > as 6 — > 0. Then, we prove that |z/ 2 1 < cp in L>i. Next, 

we prove that |u 2 (ar)| < ce^' (*) if dist(or, ri) = O(p), where | In ej = cN(S), N(S) := \ lnS\/ In | ln5\. From 
(*) Theorem 1 follows. By c, c', c, Cj various positive constants, independent of S and on T e 7a, are 
denoted. 



Step 1. Proof of the relation p{5) ^0 as 5 — > 0. Assume the contrary: 

p n := p(S n ) > c > for some sequence S n — > 0. (7) 

Let Tj n , j — 1, 2, be the corresponding sequences of the boundaries, Tj n E Due to assumption (3), one 
can select a convergent in C 2,A1 , < p. < A, subsequence, which we denote Tj n again. Thus Tj n — > Tj as 
n — > oo. From (7) it follows that (f) p(Di,D 2 ) > c > 0, where Dj is the obstacle with the boundary Tj. 
By the known continuity of the map Tj — > Aj, Tj S 7 M , it follows that Ai(a', a) — A2(a', a) = 0. 

By the uniqueness theorem [1, p. 85] it follows that Ti = r 2 . Thus, p{D\, D2) = which is a contradic- 
tion to (f). This contradiction proves that p(S) — > as S — > 0. 

Step 2. Proof of the estimate |t*2(a;)| < cp for x £ D\. It is known that ||m2||c 2 (_d^,) < c , where u 2 = 
U2(x, a) is the scattering solution corresponding to the obstacle D2- Since U2 = on T2, one has |w2(a;)| < 
(m&x xe£>i \Wu 2 \)p< cp. 



Step 3. Proof of the estimate \v(x)\ < ce dC , where v := u 2 — u\ and d := dist(x,T\). 
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From [3, p. 26, formulas (4.12), (4.17), (2.28)], one has 

\v(x)\ <e:=ceM-lN(6)}, N > *2, N{S) := i^L, 7 :=rn||>0, (8) 

0-2 > cii is an arbitrary fixed number, a,2 < \x\ < a 2 + 1 (in [3] it is assumed a 2 > aiV2, but a 2 > ai is 
sufficient). Let us derive from (8), from equation (1) for v(x), from the radiation condition for v(x), and 
from the estimate |M|c 2 (d; 2 ) < c , the estimate: 

\v(x)\ < ce d " , a; e D' 12 , c 3 p < d < c 4 p, c 3 > 0, d = dist(.T, T'J, (9) 

If (9) is proved, then Theorem 1 follows. Indeed, \v(x)\ — \v(s) + \7v ■ (x — s)\ = O(p) < ce p if d satisfies 
(9). Here we use: 1) v = u 2 — u\ = u 2 on r' l7 \u 2 \ = O(p) on r' l7 since u 2 — on T 2 , and |Vt*2| < c, 
2) \x - s\ = O(p) if dist(x,ri) = O(p), and 3) < c < |Vv| < c if d satisfies (9). The last claim follows 
from the continuity of Vtj(i), smallness of p, p(S) — > as 6 — > 0, and the fact that |Vitj|r\,- 7^ almost 
everywhere (otherwise, by the uniqueness of the solution to the Cauchy problem for (1), one concludes 
that v,j = in D'p which contradicts (2), since, by (2), \uj\ — > 1 as \x\ — » 00). Thus lnp < cp c lne, or 

(*) i^p-i-j < c/ ln(e _1 ), where p and e are small numbers, < p, e < 1, c, c' > 0, and c stands for different 

constants. It follows from (*) that p < {c/ ln^ -1 )} - ? - , where u — > as e — > 0. From the definition (8) of 
e, one gets the estimate of Theorem 1. Thus, the proof of Theorem 1 is completed as soon as (9) is proved. 

Our argument remains valid if |u| = 0(p m ) with some to, < m < 00. Such an inequality is always 
true for a solution v to elliptic equation (1) unless v = (see [26, p. 14]). 

Proof of (9). Since ||u||c 2 '''(£)i 2 ) < c > v ( x ) vanishes at infinity, and v solves (1), one can represent v(x) in 
D' 12 by the volume potential: v(x) — J Di2 g(x — y)f(y)dy, f G C^(Di 2 ), g(x) := "^j^ . The function 

|rr — 2/i = I 7 " 2 — 2r|y| cos 9+ |y| 2 ] 1//2 := R admits analytic continuation on the complex plane z = rexp{iip) to 
the sector S$ : | argz| < cj>, if z 2 — 2z\y\ cosfl+lyl 2 7^ for z in this sector. We use the branch of R for which 
ImR > 0, and ReR\j mz= o > 0. The argument of R 2 := z 2 — 2z\y \ cos 6* + \y\ 2 is defined so that it belongs 
to the interval [0, 2ir), so that the analytic continuation of g(x — y) to the sector is bounded there. It is 
crucial to have at least boundedness of the norm (f) IHIc^d^)- Indeed, (f) implies that one can extend 
v from D' 12 to Di 2 as C 1 (M 3 ) functions. This is true although the boundary dDi 2 may be nonsmooth to 
the degree which prevents using the known extension theorems (Stein's theorem, for example). The way 
to go around this difficulty is to extend u\ and U2 separately to D\ and D 2 respectively, and then take 
v = u 2 — u\ as the extension. If v € C 1 (R 3 ) satisfies the radiation condition and the Helmholtz equation, 
and is C 2 in the interior and in the exterior of D\ 2l then it is representable as a sum of the volume and 
single-layer potentials, and our argument, which uses analytic continuation, goes through. Without this 
assumption the argument is not valid and the conclusion fails, as the following example shows. 

Example 1: Let D := {x : \x\ < l,x e K 3 }, v = vi := he (1) ^ Ye(x°), where h^p{r) is the spherical 
Hankel function, Y((x°) is the normalized in L 2 (S 2 ) spherical harmonic. It is well known that h^p (r) <~ 
(^+|)r "( 2 ^"-^^~ as ^ ~~ * 00 um f° rm ly in 1 < r < 6, & < 00 is arbitrary. Therefore vi <~ r~( £+1 ) Yg (a: ) 

as I — > 00. In any annulus A := {x : 1 < 02 < r < &}, one has Ht^Hi^^) < ca 2 — > as I — > 00. On 
the other hand H'L'^ || i 2 ( s 2 ) — 1 f° r au ^- Thus, for sufficiently large I the solution to Helmholtz equation 
is as small as one wishes in the annulus A, but it is not small at the boundary 3D: for any £ its L 2 (dD) 
norm is one. The reason for the solution to fail to be small on dD is that the C 1 norm of vg is unbounded, 
as I — > 00, on dD. 
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Let us continue the proof of (9). The function v(r, x° , a), where a is the same as in (2), x° := x/r, and 
r = \x\, admits an analytic continuation to the sector S on the complex plane z, S :— {z : | arg[z — r(a;°)]| < 
<f)}, 4> > 0, r = r(x°) is the equation of the surface Ti in the spherical coordinates with the origin at the 
point O, and v(z, x°, a) is bounded in S. The angle </> is chosen so that the cone K with the vertex at r(x°), 
axis along the normal to T[ at the point r(x°), and the opening angle 2(f), belongs to D' 12 . Such a cone does 
exist because of the assumed smoothness of Tj. The analytic continuation of this type was used in [18]. 
It follows from (8) that sup r>a2 \v(r)\ < e, and sup zgS \v(z)\ < c, since Im[z 2 — 2z\y\ cos# + lyl 2 ] 1 / 2 > 
in S. From this and the classical theorem about two constants [22, p. 296], one gets \v(z)\ < ct h ^ z \ where 
h(z) — h(z, L, Q) is the harmonic measure of the set dS \ L with respect to the domain Q :— S \ L at 
the point z e Q. Here L is the ray [02, +00), dS is the union of two rays, which form the boundary of 
the sector S, and of the ray L. The proof is completed as soon as we demonstrate that h(z) ~ kd c as 
z — > r(x°) along the real axis, d := \z — r(x°)\, k = const > 0, c = const > 0. This, however, is clear: let 
r(x°) be the origin, and denote z — r(x°) by z. If one maps conformally the sector S onto the half-plane 
Rez > using the map w = z c , c' = then the ray L is mapped onto the ray L := {a?, , +00), and (see 
[22, p. 293]) h(z,L,Q) — h(z c ,L',Q'), where Q' is the image of Q under the mapping z >—> z c = w. By 

the Hopf lemma [23, p. 34], dh(0 d L ^ Q ' } > 0, h(0, V \ Q') = 0, so h(w, V , Q') - kw = kz c ' as z -> 0, and (9) 
is proved. Theorem 1 is proved. □ 



Proof of Proposition 1. It is proved in [2, p. 183] that the set {un(s, a)}v a gs 2 is complete in L 2 (T). This 
implies existence of a function is e (a,9) such that 

. f , . . , d expUO ■ s) .. 

11/ u N (s,a)is e {a,9)da — || L 2 (r) < e, (10) 

Js 2 aisl s 

where e > is arbitrarily small fixed number, N s is the exterior normal to T at the point s, and 8 € M is 
an arbitrary fixed vector. It is well known [1, p. 52], that 

-AirA{9',a) = j cxp(-i0' • s)u N (s,a)ds. (11) 

Multiply (11) by v e (a,9), integrate over S 2 and use (10), to get 

-Anlim [ A{9',a)v e (a,9)da= f cxp(-i9' ■ s) ^ff ' ^ ds. (12) 
Js 2 Jr oN s 



Note that 

dexp(i9-s) 1 f dexp[-i(6' - 9) ■ s] 



ds 



dN s 

\2 



If A 2 

= 2J y2 cx P(" iA ' x ) dx = — 2"Xd(A) 



(13) 



where the first equation is obtained with the help of Green's formula. From (12) and (13) one obtains (5). 
Proposition 1 is proved. □ 



Remark 3. In [4]- [5] attempts are made to obtain stability results for IOSP, but several errors invalidate 
the proofs in [4], [5] and [6] related to stability for IOSP. Let us point out some of the errors. Lemma 5, 
as stated in [4, p. 83], repeated as Lemma 4 in [5], claims that if a solution to a homogeneous Helmholtz 
equation in the exterior of a bounded domain D is small in the annulus R < \x\ < R + 1, \v\ <ein the 
annulus, then \v\qd < c|loge| _Cl . This is incorrect as Example 1 shows. Lemma 3 in [4] is wrong (factor 
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p 2m is forgotten in the argument). In fact, stronger results have been published earlier [17], [2], [3]. In [5] 
Lemma 2 is intended as a correction of Lemma 3 in [4] (without even mentioning [4]), but its proof is also 
wrong: the factor p 2m is not estimated. There are other mistakes in [5] (e.g., the known asymptotics of 
Hankel functions in [5, p. 538] is given incorrectly). In [6] these mistakes are repeated (p. 600). There are 
claims in [6] that: a) there is a gap in the Schiffer's proof of the uniqueness theorem for IOSP with the 
data A(a' ,a ,k) Va' G S 2 , Vfc > 0, a € S 2 is fixed [6, p. 605], b) that Theorem 6 in [8] is incorrect, and 
the proof of Lemma 5 in [8] contains a flaw [6, p. 588]. These claims are wrong, and no justifications of the 
claims are given. The remark concerning Sniffer's proof in [6, p. 605, line 1] is irrelevant (see [l,pp. 85-86]). 
It should be noted that the arguments in [4]- [5] are based on the well known estimates of Landis [9] for the 
stability of the solution to the Cauchy problem, but no references to the work of Landis are given. In [6] 
it is not mentioned that the concept of completeness of the set of products of solutions to PDE (which is 
discussed in [6] ) has been introduced and widely used for the proof of the uniqueness theorems in inverse 
problems in the works [2], [13], [19]-[21] (see also references in [2], [13]). In [24] and [25] two theorems are 
announced which contradict each other (Theorem 1 in [25] and Theorem 2 in [24]). 

Acknowledgements. The author thanks NSF for support and Prof. H.-D.Alber for useful discussions. 
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